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Holographic Implementation of the Optical Theorem 


Edwin A. Marengo, Marcus V. Bunn, and Jing Tu 


Abstract—This work provides an overview of recent advances 
in the theory and imaging applications of the optical theorem. 
The latter is a fundamental result of electromagnetics, acous- 
tics, and quantum-mechanical scattering, with many important 
applications in imaging and detection. Particular emphasis is 
given to our latest efforts aimed at developing the holographic 
implementation of the optical theorem for the measurement of 
the power extincted by a scatterer that is interrogated by an 
arbitrary probing field and is embedded in an arbitrary medium 
or scene. Several practical configurations of optical theorem 
detectors are developed in the paper. This includes a new lens- 
imaging-based optical theorem detector and different forms of 
optical theorem holograms with particular emphasis on Leith- 
Upatnieks-holography-inspired optical theorem detectors. The 
developed optical theorem holography is particularly relevant 
for optical and quantum systems in which one measures only 
field intensities, and has many potential applications including 
sensors of constitutive properties of materials, optical scanners, 
biometric readers, communication systems, and information stor- 
age devices. The optical theorem holography theory developed in 
the paper is illustrated and validated with a canonical scattering 
example in two-dimensional space. 


Index Terms—scattering theory, scattering measurements, in- 
verse scattering, holography, optical theorem, optical detection. 


I. INTRODUCTION 


The optical theorem is a basic result of wave scatter- 
ing theory, which applies to electromagnetic, acoustic, and 
quantum-mechanical scattering. In its most familiar form, 
which is applicable for plane wave excitation and free space, 
it relates the total extinction cross section of the scatterer due 
to both dissipation and scattering (into all directions) to the 
imaginary part of the forward scattering amplitude [1], [2], 
[3]. This fundamental result was first obtained by Wolfgang 
Von Sellmeier and Lord Rayleigh (John Strutt) in 1871. Lord 
Rayleigh used it in his well known investigations on the color 
and polarization of the sky. The theorem was later extended 
to quantum scattering theory by other researchers such as 
Wheeler [4] and Bohr et al. [5], and was known for some time 
as “the Bohr-Peierls-Placzek relation”. In 1952, this theorem 
was first called “a well-known theorem of optics” in a paper 
published by Rohrlich and Gluckstern [6]. It was first referred 
to as “the optical theorem” in a book published by Bethe, 
De Hoffman, and Schweber in 1955 [7]. The optical theorem 
has played a major role in the theories of electromagnetics, 
optics, acoustics, and quantum mechanics (see [8], [9], [10], 
[11], [12], [13], [14] for representative treatments), and has 


E. A. Marengo is with the Department of Electrical and Computer Engi- 
neering, Northeastern University, Boston, Massachusetts 02115, USA, email: 
emarengo @ece.neu.edu. 

M. V. Bunn is with the Department of Telecommunications, Federal 
Institute of Santa Catarina, Sao Jose, Santa Catarina, 88103-310, Brazil. 

Jing Tu is with the Department of Electrical and Computer Engineering, 
Northeastern University, Boston, Massachusetts 02115, USA. 


*This use of this work is restricted solely for academic purposes. The author of this work owns the copyright and no reproduction in any form is permitted without 


written permission by the author.* 


been at the heart of many important applications in sensing 
and imaging (see, e.g., [15], [16], [17]), and in the testing of 
computational codes ([18], [19]), among other areas. 

In recent years, both the theory and practical applications 
of the optical theorem have expanded significantly. Recent 
theoretical developments in this area include the unified 
Green’s operator theory in [20], fundamental connections to 
Green’s function extraction from field correlations [21], [22], 
generalizations to arbitrary probing fields and media [16], 
[23], [24], including extensions to nonhomogeneous media 
[25], surface waves and layered media [26], anisotropic and 
nonreciprocal media [27], open waveguides [28], transmission 
lines [29], time-varying and nonlinear scatterers [30], [31], 
active scatterers [32] and so on, as well as formulations for 
specialized sensing geometries and coordinate systems [33], 
[34], [35]. In recent years, the optical theorem has also found 
further practical applications in many fields, including lasers 
[36], metamaterials [37], sensors [38], [39], computational 
methods [40], and change detection [32], [41]. 

One of the most interesting areas that has emerged in 
recent years is the application of the optical theorem for 
detection [32], [41] and the associated generalization of the 
optical theorem to arbitrary fields and media [23], [24], [27], 
[32]. Some of these investigations [27], [32], [41] have led 
to a quite general recipe for the design of optical theorem 
sensors that measure the extincted power for arbitrary fields 
and media. This particular approach relies on the use of time 
reversal mirrors or cavities which is suitable for acoustics 
and electromagnetics. On the other hand, this method relies 
on the sensing of signal phase information, and is therefore 
not suitable for optical and quantum systems. In particular, in 
the optical regime one measures only field intensities. As a 
consequence, the corresponding optical implementation, be it 
analogically or digitally, of the required time reversal mirror, 
is less obvious. 

In this work we provide an overview of our latest efforts 
aimed at addressing this challenge. We discuss the use of 
optical holography to develop optical theorem sensors that 
measure the power extincted by a scatterer that is inter- 
rogated by an arbitrary probing field and is embedded in 
an arbitrary medium or scene that can be unknown and 
quite complex. For simplicity, we present our developments 
in the basic framework of the scalar Helmholtz equation. 
The scatterer is modeled as a nonhomogeneous and possibly 
lossy perturbation of the background which is assumed to 
be nonhomogeneous but lossless. The probing field used to 
interrogate the scatterer is general and can include, e.g., near 
field, evanescent components such as in near field imaging 
applications. Prior related work is that of Berg et al. [42] 
which illustrates an approach to measure power extinction 
using holography. However, the method in [42] holds only 
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for the typical plane wave excitation and free space scenario. 
The approach presented in this paper is more general, as it 
holds for arbitrary backgrounds and fields. Motivation for our 
research is provided by the possibility of developing novel 
sensors of constitutive properties of materials, as well as novel 
approaches for short-range optical communication systems, 
information storage devices, biometric readers, and scanners. 
Representative efforts in this area are described in [43], [44], 
[45], which address holographic fingerprinting sensors, in [46], 
[47], [48], which address hologram-based identification cards, 
and in [49], [50], [51] which discuss holographic scanning 
systems. 

It is important to point out that the proposed power sensors 
can measure both real power (due to dissipation and scattering) 
as well as the complementary reactive power associated to 
field energy storage. Importantly, the optical theorem sensor is 
based on coherent signal processing so that it exploits both am- 
plitude and phase information carried by the scattered signal. 
In this paper we show that the required coherent processing 
can be carried out analogically by means of holograms. The 
proposed approach is discussed in the classical holography 
framework ([52], ch. 9), but it can also be implemented 
via modern digital holography. It is important to point out 
that even though we focus next on the optical regime, the 
same results presented in the paper apply directly to quantum 
mechanical sensing with matter waves such as electron waves 
(where, on the other hand, instead of wave energy or power 
extinction, one deals with extinction of probability). In this 
connection, the presented optical theorem holography is also 
pertinent to change detection and tomography of fields and 
potentials using electron waves [53]. 


Il. THEORETICAL BACKGROUND 


We consider a general lossless nonhomogeneous back- 
ground medium in which propagation and scattering can 
be described in the framework of the Helmholtz operator 
[V? +k? (r)] where V? is the Laplacian and k(r) is the space- 
dependent wavenumber, which is assumed to go to the free 
space value kg = w/c = 27/X as |r| = r — œœ, where w is the 
angular oscillation frequency, c is the vacuum speed of light, 
and A is the free space wavelength. The medium response is 
characterized by the Green function G(r, r’) obeying 


[V? + k? (r)]G(r, r") = ô(r — r’) (1) 


where 6 is Dirac’s delta function, subject to the relevant 
boundary conditions. Next we focus on unbounded media 
formed by materials of finite extent so that the relevant 
boundary condition is the radiation condition. The incident 
field 7; probing the scatterer obeys 


[V? + k?(r)|pi(r) = 0 (2) 


in aregion » containing the region of interest (ROD 7 (7 C X) 
where the scatterer or medium perturbation is located. In the 
presence of the scatterer or perturbation the total field a in 4 
is given by 


Pr) = valt) + Ys(r) 


where Y, denotes the scattered field. 


(3) 


It can be shown from the main result eq. 15 of [23] that the 
average or real extincted power P. due to both dissipation and 
scattering can be expressed as the projective measurement 


P, = —kg*Im drw,(r) R(x) 
R3\r 


(4) 


where the exterior sensor or receiver mode R is such that when 
R functions as a source it generates in the ROI the complex 
conjugate (cc) probing field ;, Le., 


| dr’ R(r’)G(r,r’) = F(x) rer. (5) 
R3\r 


In addition, from the reactive optical theorem in [27], [32], it 
also follows that the extincted reactive power P, associated to 
the storage of field energy near the scatterer or perturbation is 
given by 


P, = —ky ‘Re | drw,(r)R(r) (6) 
R3\r 


where R obeys (5). 

Fortunately, there exist detectors R (indeed, an infinitude of 
them [54]) obeying the key requirement (5). This is readily 
visualized with the help of the Huygens-Fresnel principle 
which renders a singular surface source supported in the 
boundary O7 of 7 which generates the desired cc probing field 

*. This principle also provides the theoretical justification of 
the well-known time reversal method for the generation of 
the complex conjugate (time reversal) field in the background 
medium. Importantly, this approach benefits from the complex 
nature of the medium and does not require prior knowledge 
of it. In this approach the source of ; consists of a time 
reversal mirror array whose radiating elements are located far 
from the ROI at positions X,,,n = 1,--- , N and are driven by 
signals 2ikoy*ž (Xn), n = 1,--- ,N (this follows from basic 
diffraction theory and is shown, e.g., in [32], [55]). When 
this source radiates into the medium it generates a field that 
is (approximately, given the diffraction limits) proportional to 
the cc probing field * in the ROI. The theoretical principles 
and limitations of this approach are well documented in the 
literature (see [32], [55], [56] and the references therein) and 
will not be elaborated any further here. 

On the other hand, in optics, one measures only field 
intensities. Thus the optical implementation of a time-reversal- 
based optical theorem detector can be challenging. We discuss 
in the following a conceptually simple and highly practical 
holography-based approach which allows one to measure both 
P. and P, with a single optical sensor (a single-pixel camera) 
in combination with a lens and a hologram of the background 
scene. The next section outlines the basic concept for the 
simple case of a scatterer probed by an arbitrary field in free 
space, in which the detector can be based on conventional lens- 
based imaging. The following section of the paper expands 
the basic concept, by incorporating holography, making the 
technique applicable to more complex backgrounds or scenes 
that can be unknown or hidden. 
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Fig. 1. Lens-based imaging system considered in the theoretical development 
of the lens-imaging-based optical theorem detection approach. 
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Fig. 2. Optical theorem detector corresponding to a diverging probing wave 
wr in free space. 


III. LENS-IMAGING-BASED OPTICAL THEOREM 
DETECTOR 


As a first step to develop the new optical theorem hologra- 
phy, we consider the basic problem of designing an optical 
theorem detector that is based on conventional lens-based 
imaging. The developed system applies to the probing with 
arbitrary fields of a scatterer that is embedded in free space 
where the wavenumber kp = 27 /A where A is the wavelength. 
Fig. 1 shows the geometry of the system. The field used to 
probe the ROI where the scatterer may be located is generated 
by passing an incident plane wave (laser beam) through a 
space-light modulator (SLM) which can operate in either 
transmissive (as shown in the figure) or reflective modality. 
We denote the field at the object plane immediately behind 
the SLM as wo. It is assumed next that the ROI is a few 
wavelengths removed from the object plane so that the probing 
field passing by the ROI is approximately purely propagating, 
which allows us to describe propagation and backpropagation 


or time reversal in the canonical framework of the Fresnel 
approximation. Using well-known results of Fourier optics 
({52], secs. 5.3.2 and 5.3.3) we find that the field 7; at the 
image plane, where 


1 1 1 
ao (7) 
is given by 
Wr(u,v) = (u,v) ® dr(u, v) (8) 


where ® denotes convolution and A is the Fraunhofer diffrac- 
tion pattern of the lens pupil, i.e., 


h(u,v) = / / da! dy! P(Azou", Azay eT PTU Hoy’) 
— (9) 


where P is the pupil function, and w z is the ideal image field 
(under geometrical optics), i.e., 

7 1 u v 

Yr(u,v) = m” (=. a) ’ 
where the magnification M = —z2/z1 $ 

Furthermore, if the cc probing field 77 is measured and fed 

back to the imaging system (as input or “object” field at the 
original image plane), then the respective “image” field Yo at 
the original object plane is given by 


(10) 


polz, y) = h'(a, y) ® bola, y) (11) 
where 
h'(a,y) = T T dx! dy! P(Az a! , Azzy’ e? +99) 
— (12) 
and 
holz, y) = Yö lz, y). (13) 


It follows that if the effects of diffraction are neglected, then 
Wo is simply 16, the cc field at the original object plane. This 
further implies that, at least in the ideal geometrical optics 
scenario, one can reconstruct the cc probing field in the ROI. 
In practice, the diffraction effects limit the capacity to produce 
the required cc probing field exactly. On the other hand, this 
approach renders a good approximation of the required cc field 
under mild conditions. 

Assume now that a scatterer is present in the ROI. The 
scatterer is probed with the field generated by the SLM, as 
shown in Fig. 1. The objective is to estimate the total extincted 
power from intensity-only data accessible at the original image 
plane where a CCD camera can be placed, as shown in Fig. 1. 
The power budget of the corresponding scattering phenomenon 
can be measured using the optical theorem expressions (4,6). 
Note that the required receiver function R is such that when R 
is a source it radiates in the ROI the cc probing field. The lens 
analysis given above considered only the map from a field to 
another. It was found that when the field %7 is applied as input 
at the original image plane, then another field is generated at 
the original object plane which is approximately equal to the 
cc probing field Yy% . Furthermore, the corresponding field in 
the ROI is approximately equal to the respective cc probing 
field there. To apply (4,6) we need to reformulate this result in 
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terms of an equivalent source, being of particular interest one 
that is supported at the original image plane and is consistent 
with the boundary value of the field, i.e., Yï, over that plane. 
It follows from basic Rayleigh-Sommerfeld diffraction theory 
([57], ch. 2) that an equivalent source at the original image 
plane is 


R(u,v) = 2ikow7 (u,v). (14) 
Substituting this in (4,6) we obtain 
Po = -2Re | duduy; (u, v)ws(u, v) 
A 
P, = —2Im I dudvyž (u, vys (u,v) (15) 
A 


where integration is over the available sensing aperture A (e.g., 
the CCD camera area). 

In practice one employs a finite number of sensors (CCD 
detectors) or pixels as shown in Fig. 1, in which case the result 
(15) needs to be approximated as 


Pe 


K 


N 
—2KRe S > U7 (Xn)Vs(Xn) 


n=1 


P, 


l2 


N 
—2KIm )7 b7(Xn)¥s(Xn) (16) 


n=1 
where X,,,r = 1,--- , N denotes the sensor positions and K 
is a calibration constant representing the effective sensor area. 
In the absence of the scatterer, the field intensity [;(X,,) 
measured at the nth detector is (apart from noise) equal to 


In the presence of the scatterer, the output of the nth detector 
is the total intensity 


I(Xn) |w71(Xn) + s(Xn)|? 


[r(Xn)|? + Rely} (Xn): (Xn) 


where in the last statement we assume that the scatterer is 
weak so that |wW5| << |z|. If this approximation holds, then 
it follows from (16) that the value of the quantity 


N 
VOT = K» (I X 
n=1 


so that the real power P. extincted by the scatterer upon 
excitation by the given probing field can be estimated directly 
from the intensity-only data. In contrast, the measurement 
of P, via (16) is based on the same setup but requires an 
additional interferometer. 

These general results take a simple form for the particular 
case of a diverging probing field Yy = Yr produced by a point 
source or pinhole, as shown in Fig. 2. In this case, a source 
that generates the cc probing field ~* in the region where the 
scatterer is located is a point source or pinhole located along 
the optical axis at the focal point of the lens imaging system, as 
governed by the lens law. The corresponding optical theorem 
detector is the receiver counterpart of this transmitter system, 
in which the point source is substituted by a single field sensor 
or detector at the focal spot of the imaging system, as depicted 


^ 


(18) 


— (Xn)] ~ P- (19) 
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background 
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Fig. 3. Interrogation by a general probing wave of a scatterer that is embedded 
in a complex, possibly unknown background medium or scene. 


in Fig. 2. The value of the real power P. extincted by the 
scatterer is given from (19) with N = 1 by the intensity- 
only data at the single receiver. The reactive power P, can 
also be measured but only with additional interferometry. The 
basic system in Fig. 2 will remain important in the following 
developments, playing the role of the essential reference-wave 
building block of the holographic system, and allowing us to 
measure the extincted power with a single detector (a single- 
pixel camera). 


IV. OPTICAL THEOREM HOLOGRAPHY 


The main goal of the rest of the paper is to develop and 
illustrate a holographic method to measure the total extincted 
power using a single sensor (a single-pixel camera). Figure 3 
illustrates the general system under consideration, in which a 
scatterer that is embedded in a nonhomogeneous background 
or clutter is interrogated with an arbitrary probing field. Let Yp 
be the field produced by the source in the free space medium 
without the background scene. As illustrated in the figure, this 
probing field can be quite general. The field in the presence 
of the background is termed the incident field Y; and is given 
by 

Wi oa Wp 7 We 


where Ye is field scattered by the background. Finally, in the 
presence of both the background and the scatterer, the total 
field w is given by 


p = Pp + Pe + Ys = Yi + Ys 


where Y, denotes the scattered field. 

The essence of the proposed optical theorem holography 
is explained next for the special case of a probing beam 
generated by a point source, as shown in part (a) of Fig. 4(a). 
As shown in the figure, part of the probing beam can be 
relayed into a mirror for the generation of a reference wave 
Wy. As shown in part (a) of Fig. 4(a), a recording medium such 
as a photographic plate can be used to capture the intensity of 
the sum of the incident field p; plus the reference wave py. 


(20) 


(21) 
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Fig. 4. (a) Recording of the hologram of the background scene. (b) Synthesis 
of the cc probing field in the ROI using a point source or pinhole at the spot 
where the lens-based-imaging system focuses the reference wave due to the 
source and reflecting mirror in part (a). 


Under the reflective geometry assumed in the figure, Y; = we, 
so that the recording medium captures over its sensing aperture 
region A the intensity 


Iert) = |de(r) + ¥-(r)|? 
[de(r) |? +08 (r)vr (2) + ve (x) ve (2) +|vr(r)|? 
(22) 


II 


of Pe + Yr. In particular, in the following we represent the 
hologram by a transparency function ¢ over the hologram 
aperture region A that is equal to the intensity J,,, in (22) 
multiplied by a measurement constant C, i.e., 


t(r)=CIer(r) rEA. (23) 
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Fig. 5. Holographic optical theorem detector. 
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Fig. 6. Generation of the hologram using an independent reference wave 
source. 


Now, if the thus recorded hologram is excited by the cc 
reference field ~*, then it radiates back into the medium the 
cc background field = (this is the real image field of the 
hologram), as shown in part (b) of the same figure. Here we 
note that in the absence of the photographic plate, the source 
can be imaged with the lens plus detector system shown in 
Fig. 2 in which the detector position is chosen to ensure that 
the lens law is obeyed. Moreover, when the detector is replaced 
by a source, it generates in the same free space medium the 
complex conjugated version w* of ~,, as depicted in the 
figure. It follows that in the presence of the hologram, the 
cc background field ~* can be generated, approximately (up 
to the propagating field components), in the ROI, using a point 
source located at the focal point where the lens-based imaging 
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Fig. 7. Holographic optical theorem detector involving separate sources for 
the probing field and the reference field. The detector position is the (lens-law- 
defined) point where the lens focuses in free space the image of the reference 
wave source used to generate the hologram. Note that the latter source (shown 
with dashed line in the figure) is used only in the generation of the background 
hologram. This source is not used in the experiment to measure the extincted 
power by means of the optical theorem detector. 





background 


recording plate 


Fig. 8. Generation of the Leith-Upatnieks hologram for the background. The 
incidence angle of the reference plane wave is 0o. 


system images the reference wave source which in the present 
case is equal to the probing field source. In practice, it suffices 
that the thus generated field be approximately equal to the 
cc background field w= in the relevant ROI containing the 
scatterers or medium perturbations of interest. This procedure 
completes the definition of the sought-after holographic optical 
theorem detector for any scatterer in the ROI. The obtained 
holographic detector is shown in Fig. 5. Importantly, this 
holographic detector is nonadaptive or universal in the sense 
that it applies to any scatterer appearing subsequently in the 
ROI. To measure the extincted power of any scatterer in 
the ROI, the total medium including the background plus 
the scatterer in the ROI is probed using the same source 
used to excite the background during the generation of the 





background 


background hologram 


sensor 
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Fig. 9. Optical theorem detector based on Leith-Upatnieks-holography. 


hologram. The mirror or reference wave is removed, and the 
total (incident plus scattered) field is sensed using the detector 
constituted by the hologram, the lens, and the single detector. 
Moreover, this approach still holds even if the reference wave 
cannot be removed, so long as the reference beam does not 
intercept the ROI containing the scatterer. 

To simplify the exposition we have considered so far the 
special case of a point source probing the medium. This is 
a highly practical optical theorem experiment in which the 
same source can be used for probing and the generation of 
the reference wave. On the other hand, the same general 
optical theorem holography concept can be applied to arbitrary 
probing fields by adopting a separate source for the reference 
wave, as shown in Fig. 6. Of particular interest is the case of a 
point source or pinhole for the reference wave which ensures 
that the resulting optical theorem detector employs a single 
detector, as shown in Fig. 7. An important special case is 
that of plane wave excitation which corresponds to having the 
reference wave source at infinity. To simplify the discussion, 
in the remainder of this section we will focus on this particular 
special case. On the other hand, the general results presented 
next for plane reference waves remain valid for diverging 
spherical waves, as can be shown, e.g., via standard Fourier 
optics ([52], p. 119-121). 

Assume that both the background scene and the ROI are 
several wavelengths apart from the hologram aperture region 
A. Then it follows from basic Rayleigh-Sommerfeld diffrac- 
tion theory that an approximate version of the cc field %ž can 
be generated in the ROI using a source 22k9%* supported in 
the hologram aperture area A. Therefore according to (4) an 
approximate estimate P! of the extinguished power P. is given 
by 


P, ~ P, = —2Re i) dadywW(x,y)Us(x,y) 24 
A 
where here x,y denote the coordinates of point r in the area 
A of the hologram (see also the related result eq.(15)). 
If the reference wave arriving at the photographic plate is 
a plane wave w,(r) = A'etkosot traveling in the direction of 
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the unit vector sp then from (22,23) and assuming C' = 1 the 
hologram transmissivity 


t(a,y) = |e(w,y)/? + A’e0%o Ue (x, y) 


+A! e—tosor a) (ax, y) ate |.A’|? rceA. 


(25) 


For a properly designed hologram, when the background 
medium is probed with the same probing field used to generate 
the hologram, as shown in Fig. 7, the corresponding field 
intensity J, measured at the detector is given by 


Al 
Ve = SF ff acduwccenP 


where F is the lens focal length and the integration is over the 
area A of the hologram. Implicit in (26) is the isolation of the 
key real image component (the second term) in (25) so as to 
effectively reproduce the cc probing field in the ROI near the 
background when the hologram is excited by the cc reference 
field w*. This can be accomplished, e.g., via the Leith- 
Upatnieks approach ([52], sec. 9.4). In particular, referring 
to Fig. 8, if the reference wave angle 0) > sin ‘(3B)) 
where B is the bandwidth of the background, then the field 
produced by exciting the hologram with the cc reference field 
w* is approximately equal to the cc probing field ~* near 
the background (apart from a multiplicative factor which can 
be shown from (25) and the discussion leading to (24) to be 
equal to |A’|?, the fourth term in (25)). Under this condition, 
the fields due to the other three components of t in (25) (the 
first, third, and fourth terms), focus away from the ROI near 
the background (they propagate in other directions) and are 
essentially nonradiating with respect to the vicinity of the 
background. For large angles 60, the virtual image correspond- 
ing to the third term can even become evanescent. Thus as a 
radiator, the Leith-Upatnieks hologram aperture excited by the 
ce field ~* effectively radiates the cc probing field w* near 
the background, as desired. The receiver counterpart of this 
“hologram plus cc field ~*” system is the sought-after optical 
theorem detector based on Leith-Upatnieks holography and is 
shown in Fig. 9. 

Now, when there is a scatterer in the ROI the total field 
arriving at the hologram is the sum of the background field 
We. plus the scattered field ~,, and the respective intensity [.4. 
measured at the detector is given by 


(26) 


Ies =| VIe +V} (27) 
where i 
v=- [f cinean 28) 


so that if the scattered field is weak relative to the background 
field then from (26,27,28) 


|A"| 





toss Te=2 (ALY Re ff aeaniee). o 
Finally, comparing (29) with (24) we conclude that an estimate 
P” of P!, and consequently of P., can be measured from the 
intensity-only data captured by the single sensor of the optical 


theorem detector (which is located at the focal spot of the 
reference wave imaging system, see Fig. 9) via the formula 


AF 


2 
Paap =- èA) (le+s — Ic). (30) 


V. ILLUSTRATION 


We illustrate the previous developments for the system 
in two-dimensional space shown in Fig. 10. This canonical 
example models a basic reflective geometry relevant for certain 
applications such as optical scanners, biometric readers, and 
synthetic aperture radar systems, in which one is interested 
in detecting irregularities or changes over a surface (see [43], 
[44], [45], [46], [47], [48], [49], [50], [51]). In the geometry 
considered, a point source is used to probe the background 
which is modeled as a perfectly reflecting flat surface or 
mirror. As shown in the figure, the background field Ye 
consists of the superposition of the original probing field plus 
the reflected field which can be readily computed via the 
method of images. In addition, the reference field ~,. arriving 
at the hologram aperture is provided by the same point source 
used to probe the scene. In the following numerical results we 
considered a hologram aperture of length L = 400, centered 
at the origin of coordinates (x = 0,y = 0). The mirror was 
located at y = —1500X. The probing source was located at 
(c = —400\,y = —500A). The ROI was the rectangular 
region 7 = {a € R,y ER, |x + 240\| < R,y — 1500A > R} 
where FR characterizes the size of the ROI. In the simulations 
we used R = 80A. In the scattering computations we con- 
sidered a scatterer for which an exact analytical solution of 
the scattering problem is available. In particular, we used a 
perfectly reflecting semicircle of radius a < R centered at 
(x = —240\,y = —1500A). This scatterer is illustrated in 
Fig. 10, and the corresponding exact scattering computations 
for a unit-amplitude point source are provided in Appendix A. 

Figure 11 shows plots of the real part of the reflective com- 
ponent of the background field in the boundary y = —1500A 
of the ROI (solid line), and of the real part of the incident 
component of the field generated in the same boundary by the 
hologram aperture upon excitation by the cc reference wave 
(dotted line). The two plots are very similar, indicating that 
the hologram successfully generates the cc background field 
in the ROI, as desired. In the latter plot, the value of the 
field due to the hologram is multiplied by a corrective factor 
corresponding to 1/|A’|? (see the discussion following (26)), 
which is estimated numerically as the reciprocal of the average 
value of the intensity of the reference field over the hologram 
aperture. In the presence of the background mirror, the total 
field due to the hologram aperture plus cc reference wave 
source is given by the sum of this incident component plus a 
companion reflective part generated at the mirror. Figure 12 
provides plots of the imaginary part of the complex conjugate 
version of the reflective component of the background field 
in the boundary y = —1500A of the ROI (solid line), and 
of the imaginary part of the incident component of the field 
generated in the same boundary by the hologram aperture upon 
excitation by the cc reference wave (dotted line). In the latter 
plot, the field produced by the hologram is multiplied by the 


Forum for Electromagnetic Research Methods and Application Technologies (FERMAT) 


same corrective factor (1/|A’|?) involved in the companion 
(real part) plot in Fig. 11. Figures 11 and 12 convincingly 
show that the hologram, when acting as a source, effectively 
generates the cc background field in the ROI, as is required 
in order for this holographic sensor to perform the optical 
theorem operation of measuring the extincted power. 

Central to the successful holographic implementation of 
the optical theorem is the effective isolation in the ROI 
of the field produced by the real image component of the 
hologram (the second term in (25)). Figure 13 shows that this 
key requirement is achieved in this example. In particular, 
the figure provides plots of the amplitude of the incident 
fields generated in the boundary y = —1500. of the ROI 
corresponding to three components of the hologram aperture 
upon excitation by the cc reference field: The solid line plot 
corresponds to the amplitude of the real image field (second 
term in (25)). The dashed line plot corresponds to the virtual 
image field (third term in (25)). The dotted line corresponds 
to the DC or constant term associated to the intensity of the 
reference field alone (fourth term in (25)). The amplitude of 
the key real image field is significantly larger than that of the 
two other components, as desired. The field amplitude of the 
first term in (25) is not shown, but its value is comparable 
in this example to that of the virtual image. These results 
illustrate the successful isolation, in the ROI, of the key real 
image field of the hologram. This is accomplished in this 
example by putting the reference beam source away from 
the main optical axis, in analogy with the classical Leith- 
Upatnieks holography. 

The solid line in Fig. 14 illustrates the variation of the 
scattered power P, versus the scatterer radius a in the 0 — 4A 
range. The scattered power P; was evaluated using (38). 
The same figure shows the plot corresponding to the ex- 
tincted power as measured by the holographic optical theorem 
detector. Both plots are almost identical, as expected from 
Figs. 11,12,13 and the optical theorem holography theory 
presented in the paper. Figure 15 shows the percentage of 
error of the holographic optical theorem detector versus the 
scatterer radius a in the larger 0 — 160A range. The error is 
defined as error(%) = 100|P! — P!'|/ P}, where P” is defined 
in (30) while P! is defined in (38). For radius smaller than 80 
A, corresponding to a scatterer confined within the ROI 7, the 
error percentage remains below 5%. This is consistent with the 
plots in Figs. 11 and 12, which show that the cc background 
field is well approximated in this region by the transmit form 
of the optical theorem hologram. Thus this hologram obeys the 
basic requirement that allows it to function as extinction power 
sensor for scatterers residing in this ROI. On the other hand, 
as the scatterer radius increases beyond the critical radius of 
80, the sensing error grows quickly. This is expected from the 
loss of imaging capability, due to diffraction, of the finite-sized 
hologram aperture in the correspondingly larger scattering 
region, which implies that the hologram can no longer function 
(for a > 80A) as extinction power sensor. 


VI. CONCLUSION 


The optical theorem describes the energy budget of wave 
scattering phenomena. It is a basic principle that can be used 
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Fig. 10. System considered in the example. 
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Fig. 11. Plot of the real part of the reflective component of the background 
field in the boundary y = —1500A of the ROI (solid line), and of the real 
part of the incident component of the field generated in the same boundary 
by the hologram aperture upon excitation by the cc reference wave (dotted 
line). 
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Fig. 12. Plot of the imaginary part of the complex conjugate version of the 
reflective component of the background field in the boundary y = —1500A 
of the ROI (solid line), and of the imaginary part of the incident component 
of the field generated in the same boundary by the hologram aperture upon 
excitation by the cc reference wave (dotted line). 
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Fig. 13. Base-10 semilog plots of the amplitude of the incident fields 
generated in the boundary y = —1500A of the ROI by three components 
of the hologram aperture upon excitation by the cc reference field: Amplitude 
of the real image (solid line), of the virtual image (dashed line), and of the 
DC or constant term associated to the intensity of the reference field alone 
(dotted line). 
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Fig. 14. Dotted line: Plot of the scattered power Ps based on the analytical 
expression (38) versus the scatterer radius a. Solid line: Plot of the extincted 
power measured by the holographic optical theorem detector versus a. 


to design sensors of the power extincted by a scatterer that is 
immersed in a given scene or background upon excitation by a 
given probing field. In this work we have shown how to design 
optical theorem power sensors that are based on classical 
holography. The derived optical theorem holography broadens 
the perspectives on the practical implementation of the optical 
theorem. It further facilitates the development of optical theo- 
rem sensors using optical fields and quantum-mechanical mat- 
ter waves in which one measures directly only field intensities. 
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Fig. 15. Percentage of error of the holographic optical theorem detector versus 
the scatterer radius a. 


Possible applications include novel holography-based power 
sensors, detectors, scanners, communication systems, and in- 
formation storage devices based on photonics and quantum 
electronics. 

As part of the theoretical development, we detailed the lens- 
imaging-based implementation of the optical theorem in free 
space. This included the general case using arbitrary propagat- 
ing fields as well as a special case using point source excitation 
in which the optical theorem sensor can be implemented with 
a single-detector (a single-pixel camera). These results on the 
free-space lens-imaging-based implementation were applied 
later to the development of the more general optical theorem 
holography which, unlike the lens-imaging-based method, 
applies, in principle, to arbitrary scenes or backgrounds. In 
the derived optical theorem holography, information about the 
background (without the scatterer) is captured in the optical 
theorem hologram which is created by the interference of the 
background field with a known reference wave. The optical 
theorem hologram acts as the entrance coded aperture of a 
lens-based detector that is placed so as to image the reference 
wave source. The resulting hologram plus lens-based detector 
system is the sought-after holography-based optical theorem 
sensor. It is a nonadaptive sensor that measures the power 
extincted by any scatterer appearing in the neighborhood of 
the background. It performs this function well, so long as its 
transmit counterpart effectively images the cc probing field in 
the scatterer region. This was illustrated in the paper with an 
example in two-dimensional space. 


APPENDIX A 


Let (p’, ¢’) and (p, d) be the cylindrical coordinates of the 
point source and the observation point, respectively. Here the 
origin is the center of the semicircular scatterer of radius 
a. It is not hard to show using standard Green’s function 
techniques that the relevant background Green function Gp for 
the background medium bounded by the mirror, represented 
mathematically by homogeneous Dirichlet boundary condition 
for the field at the plane y = 0, is 


Galp, $; p,p) = 


i> HY (kops)Jm(kop<) sin(m@) sin(m¢’) 
m=1 
(31) 


where Jm is the Bessel function of order m, H, a ) is the Hankel 
function of the first kind of order m, p> = max(p, p'), and 
p< = min(p,p'). This Green function defines the probing 
field in the region under investigation. On the other hand, total 
Green’s function G; which corresponds to the total (incident 
plus scattered) field in the same region can be shown to be 
given by 


< Jm(koa) 1 
Gi(p,o;p\¢’) = i Jm (kop) -HO (kop<) 
2 HE (koa) 


H (kop>)sin(mø) sin(mg’). (32) 
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Thus from (31,32) the scattering Green function G's corre- 
sponding to the scattered field is given by 


Gs(p, OF p', g') = 
-i S HO (kop) HO (kop) ZE sinf(mg)sin(mg’). 
toma Hw (koa) 


(33) 


It follows from (33) and the asymptotic behavior 


z ; i 
H (kop) ~4/ “a kop > œ (34) 


that for far zone observation points 


Gs(p, Q; p, o’) oy 
. = . Im (koa) : . 1 
—iCo (i) HP (kop) Z sin(m@) sin(m¢’) 
2 HG (koa) 
(35) 
where 
Co=\ [2 cikope—in/4, (36) 
tkop 
Using (35) we can evaluate the far-zone scattered power 
P, = | pG o 
0 
ay Im(koa))? _. 
= bo AD (hop)? POO in? mg), 
> | Hm? (koa)? 
(37) 


Finally, the Bessel functions J,,(a) decay exponentially fast 
for m > x. Thus the series expansion in (37) can be effectively 
truncated, rendering the following approximate value P! of Ps: 


2 
2 [Jm(koa)] sin?(m¢’), 
|r? (koa) |? 
(38) 
where M = [koa] + p where p > 0 is a small integer. We 
used p = 5 in the computer simulations. 


M 
P, = PL = hy? XS. HQ (kop') 


m=1 
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